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TECHNICAL MEMORANDUM 


SPATIAL AND TEMPORAL ANALYSES OF ENVIRONMENTAL EFFECTS ON 
ZIZANIA PALUSTRIS AND ITS NATURAL CYCLES 


1. INTRODUCTION 


As part of a joint education and research effort funded by NASA, research studies were 
initiated involving students associated with the Ojibwe and researchers at Marshall Space Flight 
Center. Topics were chosen that satisfied the nature of the work proposed and were tractable, given 
the student’s constraints (abilities, interests, and time). One of the studies, which spanned two sum- 
mers, examined some potential environmental effects on northern wild rice in northern Wiscon- 
sin. The rice of interest is naturally occurring (‘wild’ wild rice), as opposed to cultivated wild rice 
(‘paddy’ wild rice).! 


More information on the program can be found in “Gidakiimanaaniwigamig (Our Earth 
Lodge) STEM Camp: Investigating climate change and its effect on Ojibwe lifeways,” P.I. Courtney 
Kowalczak, Fond Du Lac Tribal Community College, a proposal in response to NASA Program 
Announcement: NNH13ZHA002N-NICET-NASA Innovations in Climate Education—Tribal 
(NICE-T), funded by the NASA Office of Education, Integration, Minority University Research 
and Education Program. 


1.1 Cultural Significance of Wild Rice 


Zizania palustris (Z. palustris) (ref. 2) is known as ‘maanomim’ or the ‘good berry’ to the 
Ojibwe, who are part of the Anishinaabe culture. The Ojibwe people, termed Chippewa by the US 
government, refer to themselves as the Anishinaabe (pluralized form is Anishinaabeg). The terms 
Ojibwe, Chippewa, and Anishinaabe all have multiple spellings, and the referent peoples vary. 


The importance of wild rice to the culture of the Ojibwe is apparent even within their 
migration story. According to this, they separated from the Lenni Lenape on the east coast of the 
continental United States when Gitchi Manitou, the Great Spirit, manifested itself as a giant turtle 
in a vision to a young boy, telling him to lead his people to “where the food grows on water.” Along 
the journey, the Anishinaabe people settled along the coasts of the Great Lakes, from the western 
shores of the lower peninsula of present-day Michigan, to the eastern coast of Minnesota along 
Lake Superior. Traditionally, Z. palustris is referred to as an animate being, such as him/her, rather 
than it. Today Z. palustris is still a sacred grain to the Ojibwe. It also has economic importance to 
the tribe’s members. In 2007, the rice harvest generated more than $400 thousand in income for 
tribal members in Minnesota.? 


1.2 Study Rationale 


While northern wild rice has been extensively studied, the bulk of the studies have been 
ground or laboratory based. Remote sensing and geographic information system (GIS) technology 
has only been used to a limited extent. For example, the authors in reference 4 used Landsat data to 
map the location of wild rice in north-central Minnesota. This was a classic land cover analysis— 
1.e., where, geographically, does the targeted species occur. Minnesota Geospatial Commons has 
a similar dataset on their website for the entire state. It can be found at <https://gisdata.mn.gov/ 
dataset/env-wild-rice-lakes-rivers-wld>. The Great Lakes Indian Fish and Wildlife Commission 
(GLIFWC) has one that has a wealth of data, including wild rice covering Michigan, Wisconsin, 
and Minnesota, found at <http://maps.glifwc.org/>. Spatial analysis has also been used by authors 
such as Drewes and Silbernagel (ref. 5), who sought to understand wild rice harvesting in the con- 
text of the spatial landscape and across multiple governmental management systems. 


However, GIS technology can be used to do more than examine where something occurs. It 
can also analyze relationships between any variables that are spatially co-registered. Furthermore, 
the analyses can be performed fairly rapidly once the necessary software is minimally mastered. 
Thus, the location or extent of wild rice, if known, could be evaluated against the large number 
of spatially mapped environmental variables stored by NASA and other organizations. As many 
of the available environmental variables may be affected by changing climate, such analyses could 
contribute directly to the goals of the funded project. 


The researchers decided to test a number of environmental variables against a spatial data- 
base of wild rice, seeking potential patterns that might illuminate dependencies of wild rice on the 
climatic variables. Also, to a large degree, the research effort reported here was an effort to learn 
what the available data (which were not designed nor measured for the purposes of this study) 
could tell us if we listened. The analytical work was considered exploratory and of comparable 
importance to the educational goals. 


1.3 Biology of Zizania palustris 


The reader should note that there are two species of wild rice, Z. palustris and Z. aquatica, 
that grow in the same region. The differences between the rice are not significant here. Both are 
used by the Ojibwe in the same way, and the distinction between the species is not consistently 
maintained in nonspecialist literature. Z. palustris is an annual, subaquatic grass, active from mid- 
April through late August or early September in Wisconsin and Minnesota. For several reasons, the 
northern wild rice has been extensively studied. (See refs. 6-15 for a 10-part series of papers on the 
species’ ecology. For information on the genetics of northern wild rice, see ref. 16. Surveys of some 
of the other literature may be found in refs. 3 and 17. For a practical summary of wild rice’s ecol- 
ogy relevant to this paper, see ref. 2.) 


1.3.1 Growth 


Z. palustris has a limited range of suitable water depths, between 0.15 and 0.9 m (0.5 and 
3 ft).1” It typically grows in aquatic areas with some sort of flow, such as lakes with an inlet and 
outlet, tributaries, or small creeks and rivers. There are two critical growth stages: (1) Initial ger- 
mination (occurring in early to mid-spring), and (2) the floating leaf stage (occurring in early to 
mid-summer). In order for the seedling to germinate, a minimum 3-to-4-month dormancy period 
in near-freezing waters is needed. Germination then occurs when waters reach about 4 °C (40 °F).3 
During the floating leaf or emergent phase, Z. palustris begins exchanging gases with the atmo- 
sphere and develops buoyancy. During this stage, the rice beds are extremely susceptible to changes 
in water level and flow changes, as the rice can drown and also be uprooted by its own buoyancy 
because of its shallow root base. 


A number of things can seriously damage the rice production in a given locale. Variation in 
water level is a well-known problem. Excess inflow (precipitation or change in the upstream water- 
shed), restricted outflow (by beaver or manmade), storms, carp, other wild life, disease, insects, 
and competition from native and invasive species can all eliminate a substantial portion of the rice 
production. 


1.3.2 Boom-Bust Cycle 


It is commonly stated that northern wild rice exhibits a boom-bust cycle. Walker (ref. 18), 
the Minnesota Department of Natural Resources (ref. 3), and others—such as John Pastor and 
Tali Lee in their successful proposal to the National Science Foundation (Division of Environmen- 
tal Biology Application 0715808)—state that such cycles occur; but none of the sources examined 
offer supporting citations or substantiating data. One of the goals of this study was to substantiate 
the existence of the cycle and characterize it. 


The only reference the author has been able to find that provides data demonstrating such 
a boom-bust cycle is ref. 19. Their figure 1 shows, for the period 1970-1987, a clear boom-bust 
cycle simultaneously affecting the commercially grown wild rice harvested in both the Provinces 
of Manitoba and Ontario, while no cyclic behavior is seen in Saskatchewan. The observed cycle is 
approximately 4—S years long. It is important to note the data for the figure are production for an 
entire province. Chapter 3 discusses a province-wide drop in standard deviation of density from 
1984 to 1985 for all 20 of the study lakes in Saskatchewan. However, the subsequent 2 years exhibit 
considerable heterogeneity lake to lake. 
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Figure 1. Production of wild rice in Manitoba and price per pound in Canadian dollars. 


The Province of Manitoba has subsequently published additional statistics on wild rice 
production and price for the period 1970—2011.29 These data show that the apparent cyclic nature 
shown by reference 19 was probably not caused by a boom-bust cycle, as shown in figure 1. The 
general lack of correlation between price/pound and production indicates variations in production 
are caused by factor(s) other than the value of the harvest. 


With specific attributions, A.E. Jenks reported that multiple individuals, Ojibwe and Euro- 
pean Americans, had observed a boom-bust cycle in the wild rice.2! However, it is also clear, from 
excerpts from reference 21 (pp. 1095 and 1099) that while a 3-to-4-year boom-bust cycle existed in 
some areas, in others it did not: 


In some sections of the country the rice crop failed partly or wholly at frequent intervals. Informa- 
tion from such sources as Chief Pokagon and government farmers at Indian reservations shows that 
it so fails once in three or four years. Again, at Grass Lake, Lake County, Illinois, where there are 
1,000 acres of wild rice, it has not been known to fail in the last sixty years. 


In some sections of the country the rice crop fails partially or wholly as often as once in three or four 
years, while in other sections it has not been known to fail for long periods of time. 


Figures | and 2 of reference 22 do not strongly support a conclusion that a clear boom-bust 
cycle exists, especially not with a 3-to-4-year or 4-to-5-year cycle period. Significant annual varia- 
tion (<=40%) in northwestern Wisconsin is evident, but it does not have the amplitude of the cycles 
shown in reference 19, nor does the variation carry over into north-central Wisconsin waters. 


2. PARAMETERS OF THIS STUDY 


2.1 Data Sources 


Data for the studies reported herein are from several sources. The data pertaining to the rice 
was obtained from the GLIFWC with the assistance of their employee, Peter David. (For infor- 
mation on how the data provided by the GLIFWC are obtained, see ref. 22.) The remote sensing 
related data (air temperature and precipitation) are from the North American Land Data Assimi- 
lation System (NLDAS).23 Shape files defining the lakes and drainages of interest were obtained 


from reference 24. 


2.2 Spatial Relationships 


The study area was limited to a portion of the US termed the “Wisconsin Ceded Territo- 
ries’ (fig. 2). The Ceded Territories are lands granted from the Chippewa to the US government in 
a series of treaties signed from 1836 tol1854.2> While most of this land was granted to the federal 
government, the signatory tribes retained their hunting and gathering rights within these territories. 
This included the right to harvest the maanomin, a right retained today. 
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Figure 2. Location of water bodies included in this study and Ceded Territories 
boundary from reference 26. 


For this study, data from 40 northern Wisconsin rice-bearing lakes and riverine systems 
within the Ceded Territories were used (see figure 2 and table 1). The total east-west distance for 
the 40 water bodies is approximately 280 km (175 mi). For each of the 40 water bodies, a point 
location was obtained. Values for each water body were spatially attributed to this point. Accord- 
ing to P. David, the 40 water bodies consist of 10 small lakes (<100 acres), 10 medium lakes 
(100-300 acres), 10 large lakes (>300 acres), and 10 riverine or flowage systems (personal communi- 
cation, March 27, 2017). 


The environmental data shown in figure 2 were obtained as orthogonally gridded cells, 
which sample a continuous surface. The gridding is 1/8° of the Earth’s surface. Cell spacing in the 
study area is ~10 km on both axes. For some analyses, the gridded data were resampled to the 12 
counties in the study area. This was done by averaging all grid cells whose centroid fell within the 
county. 


Researchers recognized, prior to the work, that analyzing values reported on a per-water 
body basis against values reported on a county basis would inherently introduce a noise into the 
analyses. Furthermore, the introduced noise would have extremely complex spatial properties. 
However, the researchers felt that the magnitude of the added noise was not likely to cause a 
serious degradation of the analytical results, as the gridded surface had relatively low amplitude 
features at the scale of the counties. 


Table | is a list of water bodies included in the study. The maximum acreage for each water 
body is the largest acreage of wild rice recorded for the period 1985—2013. The maximum rice acre- 
age for some water bodies is larger than the indicated size of the water body because of differences 
in how the water body size is computed in the different databases used to assemble the table. The 
Dam column indicates whether or not the water body is a dammed drainage. A | ft dam is a dam 
that is ~1 ft high. An en-dash (—) indicates the area or status is unknown. 


Table 1. Water bodies in study. 

Max. Rice | Area 
County Name Acreage | (Acres) | Latitude | Longitude | Dam 
NF 


Barron Sweeny Pond* 40 26 | 45.392 | -91.949 
Bayfield Totogatic Lake 440 538 | 46.168 | -91.403 No 


Burnett Bashaw Lake 


Burnett Briggs Lake 


Burnett Gaslyn Lake 


Burnett Long Lake 
Burnett Mud Lake (2) 
Burnett Upper Clam Lake* 
Burnett Webb Creek* 
Douglas | Mulligan Lake 
Forest Atkins Lake 
Forest Indian/Riley Lake 
Forest Pat Shay Lake 
Forest Rat River* 

Forest Wabikon Lake 
Lincoln Alice Lake 

Oneida Fish Lake* 
Oneida Little Rice Lake (1) 
Oneida Rice Lake (3) 


Oneida Spur Lake 


Oneida Wisconsin River* 
Polk Rice Bed Creek* 
Polk Rice Lake (1) 

Polk White Ash Lakes 
Blockhouse Lake 


Billy Boy Flow* 
Blaisdell Lake 
Pacwawong Lake 


Phipps Flowage* 


Allequash Lake 
Little Rice Lake 
Manitowish River* 
Vilas Partridge Lake 
Vilas Rice Lake (4) 
Vilas West Plum Lake 
Washburn | Dilly Lake 
Washburn | Potato Lake 
Washburn | Rice Lake (2) 
Washburn | Spring Lake (1) 
Washburn | Tranus Lake 


2.3 Rice Measurements 


Two measurements and a derived value pertaining to wild rice were used in this study. The 
two measurements were rice acreage and rice density. The derived value was rice index. Rice acre- 
age is the acreage of rice on each water body. Rice density is the quantity of wild rice in a single 
acre on a scale from 0 (none present) to 5 (most abundant). P. David, who helped set up the sur- 
vey, stated the scale was intended to be approximately linear (personal communication, March 27, 
2017). A ‘S’ was designed to represent 80—100 (or more) stalks per square meter, and a ‘1’ was des- 
ignated to represent 0—20 stalks per square meter. Densities designated as 2, 3, or 4 are at 20 stalks- 
per-meter intervals. Rice index is the product of acreage multiplied by density. Therefore, rice index 
is an estimate of rice plant abundance in a single body of water. It is important to emphasize that 
rice density, acreage, and index are not equivalent to production, harvest, or yield. 


Data for the three measurements are obtained annually by visual estimation via ground and 
aerial surveys. GLIFWC has done the work for a number of years and form the longest consistent 
record of northern wild rice known to the authors. The two rice measurements are available for 
each of the 40 water bodies for each of the years 1985—2012, save for a few missing water body— 
year combinations. No information is available that numerically characterizes either the accuracy 
or precision of the measurements. As suggested by P. David (personal communication, March 27, 
2017), analyses should be sensitive to these two points. 


It is important to note that rice density is not a continuous measurement; in other words, 
a value of 2.5 is not possible. The values are binned intervals. In this respect, density is an interval 
scale.2”7 Because of the binning, it is conceptually possible that a water body that has a density of 
3 in 2 years, if more precision were available, might be 2.5 in year 1 and 3.49 in year 2. This also 
means that small numeric differences in computed values, such as average change in density, prob- 
ably have little physical significance. 


Finally, binning inherently introduces a noise compared to a continuous variable. This can 
be easily recognized if one considers the Fourier transform of a binned function versus that of 
a continuous function. The significance of this added noise depends on the details of the numerical 
operations and the magnitude of other signals in the data. 


2.4 Temporal Variables 


In addition to year, for some analyses, date ranges within a year were defined for three 
growing stages of wild rice. Germination was defined as 4/25—5/10. Floating leaf stage was defined 
as 6/10—6/25. The total growing season was defined to be 1/1—7/31. January 1 was used to capture 
snow as part of the precipitation data and for convenience, as NLDAS data are separated by year. 
No attempt was made to adjust the date ranges to compensate for either annual variations or 
geographic variation across the study area. 


2.5 Environmental Variables 


The NLDAS data were used as the source for air temperature and precipitation because the 
variables are presented in a uniform manner. The values are already surfaced, the algorithms used 
to create the data have been validated, and the data are freely, readily and publically accessible. 


Daily temperature and precipitation from the NLDAS climatological forcing data were 
obtained for an area encompassing the 12 counties containing the 40 lakes of interest. Daily tem- 
perature was obtained by averaging the 24 hourly values in NLDAS. Precipitation was the daily 
total for the cell. As previously stated, county values were obtained by averaging all grid cell cen- 
troids that fell within the county. 


2.5.1 Dammed Versus Undammed 


Water level and changes in water level are known to be very important in the biology of Z. 
palustris.3 There are a number of ways the water level in the lakes supporting wild rice can change. 
Significant precipitation in the watershed or a beaver dam, especially on a lake’s outlet are both 
possibilities. Control of water levels by manmade dams or control structures is another possibility. 
The GLIFWC provided lake damming data as GIS files for this study. Each water body was catego- 
rized as either dammed or undammed. 


2.6 Discussion of Variables 


Although only a fraction of the known parameters are considered here, the total number of 
variables in this study is still substantial. The temporal variables are year and several biologically 
relevant subsets. Spatial variables include water body size, geographic location, and various aggre- 
gation schemes, such as all water bodies in a county or subdividing water bodies based on north 
versus south or east versus west axes. Environmental variables include precipitation and tempera- 
ture with the limits of the selected temporal range. Whether a water body is subject to water level/ 
flow control is also a variable. Finally, there are the acreage and density of the rice. Thus, at the 
finest granularity, a datum has the following four dimensions: (1) time, (2) location, (3) environ- 
ment, and (4) rice. Given the constraints of time and other resources, an exhaustive investigation of 
all the variable space was not practical. What was done should be considered only a survey of the 
variable space. 


There are also an indeterminate number of uncontrolled, unmeasured, and often unidenti- 
fied, variables that affect the data. For example, there are multiple differences between each water 
body’s upstream drainage, such as the geology and agricultural usage. The species and number of 
fish can strongly affect the rice.2® Plant disease and insect infestation are active in the environment. 
Storms can destroy a water body’s rice crop. All of these, and more, add to the variation in rice 
acreage and density between water bodies and for the same water body between years. 


3. ANALYSES 


3.1 Variation in Water Body Size 


A major factor in the analysis of the data is the wide range in water body size. For example, 
the maximum rice acreage of Totogatic Lake is 22x that of Webb Creek (see table 1). This disparity 
will significantly affect any analysis where water body size is a factor. Therefore, in several of the 
analyses reported here, water body size has been normalized by one of two methods. One method 
uses the maximum reported rice acreage for the water body. The other method uses the average rice 
acreage for the water body. The former, which is the default logic in this report, has the advantage 
that all normalized values are between 0 and 1. Its disadvantage is that the annual values inherently 
tend to be in a limited portion of the available dynamic range. The latter method, using the average 
for normalization, has the advantage that the dynamic range is numerically wider, but the dynamic 
range varies between water bodies. The median rice acreage cannot be used for normalization, as 
some water bodies have 0 acreage for more than half of the years. 


3.2 Annual Variation of Rice Acreage and Density Among Water Bodies 


The range of normalized acreage for most years is 0-1, meaning that in any given year, at 
least one of the water bodies is at its maximum reported rice acreage, and at least one water body 
had no rice acreage. Only 9 years have a maximum value below 1, and all 9 of those years have 
maximum values 20.83. Two years have a minimum greater than 0. 


Standardized deviation of the normalized rice acreage was obtained, see figure 3. The lim- 
ited range of standard deviation from year to year strongly demonstrates that the sampled 40 water 
bodies have similar variation between water bodies within a year. It is noteworthy that the standard 
deviation has not changed noticeably while the normalized average acreage has dropped by at least 
half. 
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Figure 3. Average over all water bodies of annual normalized rice acreage 
and standard deviation. 


The range of density is similar in behavior to the range of acreage, uniformly 0—5 in almost 
all years. The standard deviation of density (fig. 4), is also similar in nature to the standard devia- 
tion of acreage—fairly consistent irrespective of a long-term drop in average density, which is shown 
in figure 5. The variation of rice acreage and rice density between water bodies in a single year is 
very large. This is a major signal in the dataset. 


Standard Deviation of Density 
ipo} 
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Year 


Figure 4. Standard deviation of rice density for all water bodies within a year. 
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Figure 5. Average rice acreage and density of the water bodies. 


3.3. Interannual Variation of Rice Acreage and Density Within Water Bodies 
The interannual variation in rice acreage (fig. 6) is computed by the following algorithm: 
(1) The water body’s average rice acreage over the time period 1985-2012 is computed. 


(2) For each water body, the annual acreage is divided by the water body’s average rice acre- 
age. This removes the effect of large versus small water bodies. 


(3) The difference between the normalized acreages for sequential years is then computed. 
(4) For each water body, the average of all the annual differences is then obtained. 


There is large interannual variation in the rice acreage for the water bodies. This may be the 
most striking pattern in the data. Eighty-five percent of the water bodies in the study have an aver- 
age interannual variation between 30% and 70%. In other words, it is common for rice acreage to 
change by 50% year to year. Sweeny Pond typically varies more year to year than the average of its 
rice acreage. The annual acreage of Totogatic Lake (marked ‘Largest’ in fig. 7) varies year to year 
by 89%. 
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Figure 6. Frequency of the average variation in rice acreage for a water body. 
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Figure 7. Average normalized rice acreage for water body quartiles, which are defined based 
on maximum rice acreage for each water body. 


There is no obvious pattern to the interannual variation of acreage based on lake versus 
flowage or the size of the water body. This intense interannual variation obscures most other varia- 
tions, and it is not explained by the variables examined in this study. 
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The existence of a boom-bust cycle can be readily checked using the interannual variation 
data and the following logic: Compare the acreage in a water body in year | to the acreage for the 
same water body in year 2. If the acreage increases in the second year compared to the first year, 

a flag value of 1 is recorded. If instead the acreage decreases, a flag value of —1 is recorded. If no 
change between the years is observed, a flag of 0 is reported. When this operation is done for the 
40 water bodies for all the years, the sum of the flag values is 1,094 (12 year-water body data points 
are missing), is +20. 


Compare this to what would be expected if a 4-year boom-bust cycle existed. For a single 
water body, for each cycle, three 1s would be recorded and one —1, with a sum for the cycle of +2. 
A 3-year cycle would have a sum of +1. Using a 4-year cycle, if it existed consistently, over all the 
water bodies for all the years, the global sum would be 7 (number of possible 4 year cycles) x 40 
(number of water bodies) x 2 (value of a cycle) = 560. If interannual variation in acreage was ran- 
dom, the value would tend to 0. These data strongly support a hypothesis that interannual varia- 
tion is effectively random, rather than a hypothesis that a multiyear boom-bust cycle exists. 


There is also interannual variation in rice density. A check for a possible boom-bust cycle 
in density can be done in exactly the same way as for acreage. When this is done, the sum is +31, 
which is essentially what would be expected from random interannual variation. The magnitude of 
the average interannual density variation is small. Thirty-eight of the water bodies have an aver- 
age variation less than 0.18 units. The other two water bodies have average variations of 0.35 and 
0.25 units. 


A check was made to see if acreage and density vary in a correlated manner. This was 
done by comparing the corresponding flag values for each water body for each year. If both flags 
were +1 or both were —1, a second flag of +1 was recorded. If the acreage and density flags had dif- 
ferent signs, a second flag value of —1 was recorded. If the flag for either density or acreage was 0, 
a second flag of 0 was recorded. With this logic, the second flags were then summed. It was found 
that out of 1,093 water body—years, rice density and rice acreage both increased, or both decreased 
486 times. Either density or acreage stayed the same 377 times. The density changed with an oppo- 
site sense of the change in acreage only 230 times. This is strong support for the conclusion that 
density and acreage do co-vary, but the relationship is not extremely strong. 


3.4 Multiyear Variation in Rice Acreage and Density Within Water Bodies 


In addition to the interannual variation, there are variations in rice acreage and density that 
occur on longer time spans. Figure 5 shows the average acreage of wild rice in the surveyed water 
bodies versus year. Nineteen of the years have data for all 40 water bodies. Seven of the years have 
data for 39 water bodies. The remaining 3 years have 38, 37, and 36 water bodies. 


There is no evidence of a boom-bust cycle when the data are examined this way. However, 
there are several very important observations that can be drawn from the two curves: Over the 
almost three decades of measurement, the expected acreage has dropped by approximately 50%. 
Density also clearly shows a long-term decrease, though the magnitude of total decrease in density 
is not as large as that seen in the acreage data. Another observation is there can be strong (=50%) 
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variations in acreage over 4-year periods. Although there are variations in density over time peri- 
ods shorter than a decade, density clearly does not vary in the same way as acreage or in synchrony 
with acreage. 


Points | and 2 together suggest that at least two phenomena can substantially affect acreage. 
One of the phenomena might or might not also affect density. Whatever the phenomena are, they 
can operate over hundreds of kilometers. 


Approximating the large swings in average acreage over multiple years, e.g. 1987-1990 and 
1997-2000, by the summation of a random variable acting on the 40 water bodies is possible; 
however, this possibility was not quantitatively investigated. Even if this were regarding this data, it 
certainly would not seem to apply to the rice harvest in Manitoba, which shows even larger, sharper 
swings (fig. 1). 


Domination of the average acreage by a few large water bodies is also conceptually possible. 
For example, at its maximum coverage, Totogatic Lake in Bayfield County contains 1/8 of the total 
maximum rice acreage (see table 1). This possibility can be checked by normalizing the annual acre- 
age in each water body by the maximum acreage observed for that water body. This was done, and 
the average normalized acreage plotted versus year (fig. 3). Similar patterns are present in figure 3 
as in figure 5, though they are reduced in amplitude because of the method of normalization (using 
the maximum rice acreage for each water body, thus values are constrained to the range 0-1). This 
suggests that a few large water bodies are not driving the observed patterns. As the dynamic range 
of density is 0-5 in steps of 1, normalization of density was not needed. 


The phenomena can possibly acting on the average water body acts differently as a func- 
tion of water body size. This was checked by grouping the water bodies into quartiles and plotting 
the average normalized acreage for each quartile (see fig. 7). There are four water bodies in the 
75-to-100 rice acreage quartile, five in the 50-to-75 quartile, 10 in the 25-to-50 quartile, and 21 in 
the 0-to-25 quartile. Also shown is the normalized rice acreage for Totogatic Lake, which has the 
largest maximum rice acreage of the 40 water bodies. Note the swings in acreage that one water 
body exhibits. There is no clear difference in normalized average rice acreage that is related to the 
maximum rice acreage. 


3.5 Environmental Variables 


For this work, daily average temperature and daily precipitation for each NLDAS grid cell 
were obtained for the entire year. Annual mean precipitation and annual mean temperature were 
obtained for each county by averaging all the daily precipitation or temperature values for all grid 
cells whose centroid fell within the county. This process obtains an approximation of a more rigor- 
ous area weighted average while simplifying processing details. 


For each year, seasonal mean temperatures for each of the three seasons—germination, 


floating leaf, and growing—were obtained for each county by averaging all the daily temperature 
values during the season for all grid cells whose centroid fell within the county. 
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The three seasonal precipitation values data were handled differently. As for temperature, 
a daily precipitation value for the county was obtained by averaging all the daily precipitation value 
for all grid cells whose centroid fell within the county. Then, for each year and each county, the 
daily precipitation for each of the three seasons were summed over the season’s time period. This 
gave a season total precipitation by county for each year. 


In total, there are eight values for each county for each year—annual mean precipitation 
and temperatures, season mean temperatures for three seasons, and season precipitation summa- 
tions for three seasons. Values are obtained by the methods explained in sections 2.5 and 3.6. To 
permit plotting on a single axis, the daily average values were multiplied by 271, and the growing 
season values were multiplied by 0.128. 


As stated previously, the variations in rice acreage within a year and between years are 
both very large. In contrast, variation of precipitation and temperature across the 12 counties in 
a single year is much smaller than the variation for a single county across the 28 years, 1985-2012 
(see table 2). Stated another way, the variation between years for a county is much greater than the 
variation within a year between counties. Any analysis using values by county would have substan- 
tially more noise than an analysis using values by year. Therefore, analyses were done on a yearly 
basis. 


Table 2. Comparison of average standard deviation by county versus by year. 


Floating 
Germination Leaf | Growing 
Season Season | Season | Annual 


Standard Deviation in Precipitation 
For 12 counties 23.6 31.1 99.4 0.37 


For 28 years 10.1 13.5 37.8 0.15 
Standard Deviation in Temperature 

For 12 counties 44 3.6 2.7 2.06 
For 28 years 13 1:5 0.8 0.65 


Figures 8 and 9 show the four precipitation and four temperature variables examined in this 
study versus year. None of the eight plots show clear trends or patterns with time. The fact that 
there is no spatial component to these plots is emphasized. 
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Figure 8. Annual variation and long-term trend in precipitation for the daily average, 
germination season, floating leaf season, and growing season. 
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Figure 9. Annual and long-term variation in temperature for the daily average, 
germination season, floating leaf season, and growing season. 
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3.6 Environmental Variables Versus Rice Variables 


The relationships between temperature or precipitation versus rice acreage, density, or index 
were investigated by examining the annual average of each of the three rice measures across all 
water bodies versus the corresponding yearly average or seasonal summation precipitation or the 
average temperature. This made a total of 24 graphs. Figures 10 and 11 are typical of many of the 
graphs that do show a relationship. Again, emphasizing these graphs show no spatial relationships; 
but they show time- environment-rice is important. 


Table 3 shows which of the 24 graphs visually exhibit relationships. None of the graphs 
show notably stronger relationships than those exhibited in figures 10 and 11. When a relationship 
is seen with rice acreage, acreage is increasing with decreasing precipitation or decreasing tempera- 
ture. When a relationship is seen with rice density, density is increasing with decreasing temperature. 


Average Acreage 


0 200 400 600 800 1,000 1,200 1,400 1,600 
Summed Precipitation (mm) 6/10-6/25 


Figure 10. Average rice acreage over all 40 water bodies versus the sum of precipitation 
during the floating leaf period. The dashed line is the least squares linear fit 
to the data. R? is 0.22 for the line. 
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Average Rice Acreage 


Average Annual Temperature (°F) 


Figure 11. Average annual temperature for all counties versus average rice density 
for all water bodies. R2 = 0.24. 


Table 3. The presence/absence of an apparent relationship between precipitation 
and temperature for various time periods versus rice acreage, density, 


and index. 
Time 
Period Acreage | Density 


Precipitation | Annual 


Growing 


Germination 


Floating Leaf 
Temperature | Annual 


Growing 


Germination 


Floating Leaf 


Density shows only two occurrences of a possible relationship. It is certain that no strong 
relationship exists. It may also be true that weaker relationships may be obscured by the binned 
nature of the variable with its limited dynamic range. 


3.7 Effects of Damming 
Preliminary analysis appears to show that grouping dammed and undammed water bodies 


produce no significant difference results. Refinements to the methodology should be considered 
before completely ruling out manmade water level controls as a factor on these waters. 
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4. DISCUSSION 


4.1 This Study 


The choice of variables for this study was largely dictated by data availability and time 
constraints, not by the biology of Z. palustris. For example, biologically, it is highly probable that 
during germination, it is water temperature—not air temperature—that is relevant; but the latter is 
available, and the former could only be approximated with modeling that probably has never been 
done. 


The existence of other more meaning variables that have been measured or can be estimated 
is entirely reasonable. For example, wild rice self-seeds by shattering. The timing and nature of 
wind could play a significant role in the overall process, and some measurements of wind can be 
recovered from the meteorological databases. It is also entirely possible that local storms could 
play a significant role in the variation between water bodies in a single year, and in the interannual 
variation for a single water body. Storms can damage or destroy the plants during the floating leaf 
stage. Conceptually, storms later in the growing season could cause lodging, which is a well-known 
problem with many seed-bearing grasses, such as wheat. In either case, storm damage could greatly 
reduce the seed supply for the subsequent year. Also, figure 11 and an anecdotal observation from 
P. David (written communication from April 26, 2017) suggest that colder winter temperatures— 
which may be obtained from the NLDAS data—could be important, as they might cause better 
seed germination. 


The two biggest signals in the rice data are clearly the variation among water bodies within 
a year and variation for each water body between years. Metaphorically, these are the elephants 
in the room. This study has not identified any evidence as to the nature of the processes that are 
responsible for these variations. Within the limits of this study, it is possible that a single process or 
a collection of processes causes both. We also emphasize that the strength of the two big signals in 
the rice data make detection of the weaker signals much more difficult. In fact, it would have been 
impossible to demonstrate the existence of the weaker signals without the large, long-term dataset. 


The lack of a strong relationship between interannual changes in rice acreage and rice den- 
sity 1s striking. It means that very large changes in area covered by rice occur without correspond- 
ing changes in density. As the coverage of wild rice for a water body generally returns at some 
subsequent year, it is clear the interannual change is not due to degradation/destruction of the rice 
bed. This is another major signal in the data without known causative agent(s). 


The subtler signals in the data are also unexplained. These include the long-term decrease in 


acreage and density in the study area and the relationships between some combinations of tempo- 
ral-precipitation and temporal-temperature versus rice acreage or rice density. 
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4.2 Other Datasets 


For a single water body, the changes in rice acreage year to year raises questions about the 
spatial patterns of the changes. We point out that knowledge of the changing spatial patterns on 
an interannual basis could be useful in building or restricting hypotheses about the cause(s) of 
the interannual variations. Vogt refers to a library of imagery for approximately 80 water bodies 
in Minnesota covering the time period 2007—2016.?? P. David has published multiple aerial pho- 
tographs of rice-bearing waters and has shared some impressive time sequence images with the 
author (personal communication 4/26/2017). There may be other archives of aerial photographs 
that would be useful. 


Water depth is a major ecological variable for Z. palustris. While water depth can be mea- 
sured multiple ways from a boat, surveying a significant number of water bodies that way becomes 
a sizeable task. For the type of lake of interest for wild rice, it 1s very likely this variable can be 
accurately recovered from multispectral imagery obtained by various satellite-borne instruments. 
Clear water has moderate transmissivity in bluer wavelengths and decreasing transmissivity into the 
near-infrared. Analysis of imagery acquired when there is no vegetation on the water body surface 
can use this change in transmissivity to estimate depth of the water.9 


There are a number of meteorological variables that can be obtained in surfaced form or 
can be readily surfaced. Lightning is an excellent proxy for thunderstorm activity. Lightning is mea- 
sured from space by such instruments as the Lightning Imaging Sensor on the International Space 
Station, and the Geostationary Lightning Mapper on GOES-16. The National Oceanographic 
and Atmospheric Administration archives all of the ground-based meteorological data for first- 
order and co-op stations. Measurements include temperature, humidity, winds, precipitation, cloud 
cover, and visibility. All of these are hourly data. In addition to temperature and precipitation, as 
used here, NLDAS also has wind speed, wind direction, and specific humidity (from which rela- 
tive humidity can be computed). NLDAS land surface model outputs could also be helpful; these 
include surface runoff, total evapotranspiration, soil moisture at four vertical depths, soil tempera- 
ture at four vertical depths, vegetation fraction, and snow cover fraction. 


4.3 Other Rice Data for Future Work 


For future work, study of other wild rice datasets could be valuable inclusions. 
One such dataset is in reference 29. Data and analyses for 10 lakes measured over the time period 
1998-2016 is presented. The publication includes information on multiple variables not readily 
available or nonexistent for the GLIFWC data used here. The variables include water depth, water 
temperature, and water quality. The rice parameters measured are density, average stalks per square 
meter, and biomass, in grams per square meter. One interesting observation in this work shows 
a decrease in rice biomass over the survey period but not a decrease in acreage. 


A second source of data is rice production in Manitoba. The production data (fig. 1) for 
Manitoba do not show a long-term decline. It could be there is no long-term decline in Manitoba; 
equally, some change—such as increasing harvesting activity—is masking a decline. The informa- 
tion included here does not illuminate the situation. The Manitoba data does show wild annual 
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swings in production. This is Province-wide data, which suggests that whatever is causing such 
intense changes in production operates on scales of many hundreds of kilometers. There are vari- 
ous management differences between the Manitoba rice waters and Wisconsin waters, of which any 
joint analysis would need to be aware. 


Genetically, the paddy wild rice is essentially identical to the wild wild rice, though it has 
a slightly smaller range of characteristics due to selection for characteristics such as resistance to 
shattering. It would be prudent to investigate the observations and production records of commer- 
cial wild rice producers in the region. By creating fields, which are inherently consistent relative to 
the open, natural system studied here, the farmers may eliminate or control variables which cannot 
be isolated in a study such as this. It is possible both growers and researchers and the Ojibwe who 
depend on the wild rice harvest could all benefit from such a collaboration. 


4.4 Statistical Significance 


Note that this report has purposely avoided the use of parametric statistical analysis to 
make assertions about significance. This is for several reasons. First, the nature of the underly- 
ing distributions is unknown. Secondly, application of the central limit theorem in this case is at 
best problematic, as it applies to independent, random variables. It seems highly likely that there 
are multiple dependent, nonrandom variables involved in this system. Whether any or all of the 
system can be approximated by random variables has not been established. Thirdly, one of the 
major parameters, rice density, is a binned, low-resolution measure. While the numerical operations 
required to compute statistical significance can be performed on such data, the meaningfulness of 
the result may not be as great as expected. Finally, we repeat the axiom ‘Correlation is not proof 
of causation.’ This study has identified a large number of relationships. It has not identified any 
causative or controlling forces. Asserting significance inherently tends to obscures this. 


4.5 Summary of Observations 


To assist the development of future research hypotheses, the major observations of this 
study are summarized here. An understanding of the ecology of Z. palustris should help explain 
the following observations: 


(1) The annual variation of rice acreage and rice density among water bodies is very large. This is 
a major signal in the dataset. 


(2) The consistency of the annual range of rice acreage and density, as well as the limited range 
of standard deviation from year to year, strongly demonstrates that the sampled 40 water bodies 
have similar variation between water bodies within a year. 


(3) It is noteworthy that the standard deviations of rice acreage and rice density have not changed 


noticeably, while the normalized average acreage has dropped by at least half, and the average den- 
sity has also dropped substantially. 
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(4) There is large interannual variation in the rice acreage for the water bodies. This may be the 
most striking pattern in the data. 


(5) There is no obvious pattern to the interannual variation of acreage based on lake versus flow- 
age or the size of the water body. 


(6) Rice acreage can change radically without comparable changes in rice density. 


(7) It is clear these data strongly support a hypothesis that interannual variation is effectively ran- 
dom, rather than a hypothesis that a multiyear boom-bust cycle exists. 


(8) The data strongly support the conclusion that density and acreage do co-vary, but the rela- 
tionship is not extremely strong. 


(9) Over the almost three decades of GLIFWC measurements, the expected rice acreage has 
dropped by approximately 50%. Density also clearly shows a long-term decrease, though the mag- 
nitude of total decrease in density is not as large as that seen in the acreage data. 


(10) There can be strong (=50%), noncyclic variations in average acreage over 4-year periods. 


(11) Although there are variations in density over time periods shorter than a decade, density 
clearly does not vary in the same way as acreage or in synchrony with acreage. 


(12) Whatever phenomena control the long-term trends in average rice acreage and density oper- 
ates over hundreds of kilometers and many years. 


(13) A few large water bodies are not driving the observed patterns. 


(14) The maximum rice acreage of a water body does not appear to have a significant relationship 
to the water body’s normalized average rice acreage over time. 


(15) When a relationship is seen between rice acreage and air temperature or precipitation, it 1s 
such that acreage increases with decreasing precipitation and decreasing temperature. 
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APPENDIX A—RICE ACREAGE 


Data in tables 4-6 are courtesy of Peter David, Great Lakes Indian Fish and Wildlife Com- 
mission. As the data are visually estimated, there is the possibility of human variation in the mea- 
surements. Mr. David comments in a written communication to D. Rickman (April 26, 2017) that 
the data for 1985 and 1986 were obtained by his predecessor and that work may have influenced his 
own work in the early years. Removal of this data from the analyses would not substantially affect 
the results of this study. 


In 2007 Upper Clam Lake developed a problem with carp, possibly due to a loss of blue gill 


(Havranik, 2012). Intensive restoration efforts have been undertaken since ~=2011. Mulligan Lake 
has had a beaver dam for several years, which has destroyed the rice habitat. 
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Table 4. Rice acreage for each lake for years 1985-1995. 


County Water Body 


Barron Sweeny Pond 
Bayfield Totogatic Lake 
Burnett Bashaw Lake 


Burnett Briggs Lake 


Burnett Gaslyn Lake 


Burnett Long Lake 
Burnett Mud Lake (2) 
Burnett Upper Clam Lake 
Burnett Webb Creek 


+~+|O}COlNI| OD] a] &{]w]rm]— 


Douglas Mulligan Lake 


x 


Forest Atkins Lake 
Forest Indian/Riley Lake 
Forest Pat Shay Lake 
Forest Rat River 

Forest Wabikon Lake 
Lincoln Alice Lake 

Oneida Fish Lake 

Oneida Little Rice Lake (1) 
Oneida Rice Lake (3) 
Oneida Spur Lake 


x 


Oneida Wisconsin River 
Polk Rice Bed Creek 
Polk Rice Lake (1) 
Polk White Ash Lakes 
Blockhouse Lake 


Billy Boy Flow 
Blaisdell Lake 
Pacwawong Lake 


Phipps Flowage 


Allequash Lake 
Little Rice Lake 
Manitowish River 
Vilas Partridge Lake 
Vilas Rice Lake (4) 
Vilas West Plum Lake 
Washburn | Dilly Lake 
Washburn | Potato Lake 
Washburn | Rice Lake (2) 
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Table 5. Rice acreage for each lake for years 1996-2007. 


County Water Body 1996 | 1997 | 1998 | 1999 ) 2000 | 2001 | 2002 | 2003 | 2004 | 2005 | 2006 | 2007 


Barron Sweeny Pond 7 15 8 3 5 3 5 | 20 1 11 0 1 
Bayfield Totogatic Lake 15 | 440 | 135 | 95 | 51 65 | 18 | 120 | 135 | 350 | 108 | 215 
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Forest Wabikon Lake 
Lincoln Alice Lake 
Oneida Fish Lake 
Oneida Little Rice Lake (1) 
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Oneida Spur Lake 
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Little Rice Lake 
Manitowish River 
Vilas Partridge Lake 
Vilas Rice Lake (4) 
Vilas West Plum Lake 
Washburn | Dilly Lake 
Washburn | Potato Lake 
Washburn | Rice Lake (2) 
Washburn | Spring Lake (1) 
Washburn | Tranus Lake 


Table 6. Rice acreage for each lake for years 2008-2013. 


County Water Body 2008 | 2009 | 2010 | 2011 | 2012 | 2013 


Barron Sweeny Pond 1 8 3 11 3 0 
Bayfield Totogatic Lake 54 | 180 | 81 | 110 | 35 | 58 
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Manitowish River 
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Washburn | Rice Lake (2) 
Washburn | Spring Lake (1) 
Washburn | Tranus Lake 


APPENDIX B—RICE DENSITY 


Data in tables 7-9 are courtesy of Peter David, Great Lakes Indian Fish and Wildlife Com- 
mission. For comments, see Appendix A. 
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Table 7. Rice density for each lake for years 1985-1995. 
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Table 8. Rice density for each lake for years 1996-2007. 


County Water Body 1996 | 1997 | 1998 | 1999 2004 | 2005 | 2006 | 2007 
0 4 
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Table 9. Rice density for each lake for years 2007-2013. 
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